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ABSTRACT 

We present a photometric study of star clusters in the nearby starburst galaxy 
M82 based on the UBVI, YJ and H band Hubble Space Telescope images. We find 
1105 star clusters with V < 23 mag. 1070 of them are located in the disk region, 
while 35 star clusters are in the halo region. The star clusters in the disk are 
composed of a dominant blue population with a color peak at (B — V)q ~ 0.45, 
and a weaker red population. The luminosity function of the disk clusters shows 
a power-law distribution with a power law index a = —2.04 ± 0.03, and the scale 
height of their distribution is h z = 9". 64 ± 0".40 (164 ± 7 pc), similar to that 
of the stellar thin disk of M82. We have derived the ages of ~630 star clusters 
using the spectral energy distribution fit method by comparing UBVI(YJ)H band 
photometric data with the simple stellar population models. The age distribution 
of the disk clusters shows that the most dominant cluster population has ages 
ranging from 100 Myr to 1 Gyr, with a peak at about 500 Myr. This suggests 
that M82 has undergone a disk-wide star formation about 500 Myr ago, probably 
through the interaction with M81. The brightest star clusters in the nuclear 
region are much brighter than those in other regions, indicating that more massive 
star clusters are formed in the denser environments. On the other hand, the 
colors of the halo clusters are similar to those of globular clusters in the Milky 
Way Galaxy, and their ages are estimated to be older than 1 Gyr. These are 
probably genuine old globular clusters in M82. 

Subject headings: Galaxies: individual (M82) — Galaxies: irregular — Galaxies: 
starburst — Galaxies: star clusters: general 
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Introduction 



Nearby starburst galaxies have been of great importance in terms that they provide 
a detailed view of the starburst process in galaxies. Especially, thanks to the high spatial 
resolution images obtained with the Hubble Space Telescope (HST), it became possible to 
study the properties of the star clusters and to use them for investigating the star formation 



or star-forming galaxies (e.g. 


Whitmore et al. 


1999 


Larsen 


2002 




Scheepmaker et al. 


2007; 


Park et al. 


2009; 


Hwang & Lee 


2008 


Cantiello et al. 


2009; 


Hwang & Led 


2010 


Pellerin et al. 


2010; 


Whitmore et al. 


2010; 


Chandar et al. 


2010[ 


Larsen et al. 


2011 


). It is found that young 



star clusters in these galaxies sh ow a power law luminosity funct i on (N cx My) with s imilar 
power index values (a ~ —2) ( IBastian. Covey, fc Meyer 1 12010| ; iHwang &; Led 120081 ) , and 
that some of them are more massive (1O 5 M < M < 10 6 M Q ) than normal young clusters 
(M < 10 4 M o ) so that they are called super star clusters (SSCs). 

Existence of an old halo in these galaxies also has been known. Several studies in- 
vestigated old star clusters in star-forming or starburst galaxies. A significant population 
of old star clusters is found even in late-t ype galaxies with high star formation rate (e.g. 
Chandar et al.ll2004J ) . IHwang fc Led ( 120081 ) found several tens of old globular cluster candi- 
dates in an interacting and highly star- forming spiral galaxy M51, and estimated its globular 
cluster specific frequenc y to be as low as ~ 0.25, which is even lower than that of the 
Milky Way (Sn ~ 0.5). IWhitmore et al.l ( 120101 ) reported the presence of some old clusters 



in the Antennae galaxy with high star formation activity. 



M82 is one o f the nearest starburst galax i es, and its star format ion rate is as high 



as ~ 



1OM yr" 1 f fo'Connell fc Manganol Il978t iBell fc Kennicuttl l200lh . one of the high- 



est values in the nearby universe (see Mayya fc Carrascol ( 20091 ) for an extensive review). 



ral arms in the near-infrared (N IR) images led it to be classified as SBc ( Telesco et al. 



(Telesco et al. 


1991 




Greco, Martini. & Thompson 



20121 ). This galaxy belongs to the M81 Group and it is embedded in a huge neutral hydroge n 



cloud connecting all major member galaxies ( lYun. Ho. fc Lolll994j ; IChynoweth et al.l 120081 ). 



The observation of group-wide neutral hydrogen structure provided a strong evidence that 
the interaction between M81 and M82 might have induced the starburst phenomenon in M82 
feieke et al.ll98pl : lYun. Ho. fcTolll994i lYuJl999l : iForster Schreiber. Genzel. Lutz. fc Sternberg 



2003 



There have been many studies on the star clusters in M82 since lO'Connell fc Mangano 



( 119781 ) reported the existence of several knots and semi-stellar objects in the central region 
of M82. Some of them were classified as SSCs based on their high luminosity (—13.2 < 
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M B < -9.6) (IQ'Connell et all Il995h . Later Melo et all fl2005[ ) detected 197 SSCs with 
17 < F555W < 21 in the nuclear region of M82 using the data taken with H ST /Wide Field 
Planetary Camera 2 (WFPC2), and concluded that m ost nuclear star clusters in M82 are 
young (age < 25 Myr) and compact (r = 5.7 ± 1.4 pc). iBarker. de Grijs. fc Cervind (120081 ) 
reported an enhanced star cluster formation at about 100 Myr ago in the inner few hundred 
parsecs of M82 based on the age distribution of about 150 star clusters. 

Other than t he central or nuclear region, a post-star burst region B has been a target 
of many studies. Ide Grijs. O'Connell. fc Gallagher! (120011 ) detected 113 star clusters in the 
region B using the data taken with ifST/WFPC2 and Near Infrared Camera and Multi- 
Object Spectrometer (NICMOS) and concluded that the last tidal encounter between M81 
and M82 occurred at about 600 M y r ago based on the age distribution of the star clusters. 
Later de Grijs. Bastian. Sz Lamersl ( 2003 ) revis ited the same data set and argued that the 
encounter epoch is about 1 Gyr ago. However, ISmith et al.l (120071 ) found 35 star clusters in 
the region B using the image data obtained with i/ST/Advanced Camera for Surveys (ACS) 
and concluded that M81 and M82 might have encountered at about 150 Myr ago based on 



the ag es of the detected star clusters. This result is supported by Konstantopoulos et al. 



(120081 ) who showed that the spec troscopic age s of se ven star clusters in the region B are in 
agreement with those reported in ISmith et al.l (120071 ). 



A study of the star clusters over the entire disk of M82 was given by iMayya et al. 



(120081 ) who found 653 star clusters in the disk of M82 using the image data obtained with 
H ST j ACS. However, they assumed two fixed values for the ages of star clusters, 8 Myr 
and 100 Myr, instead of deriving their ages, to study cluster mass distributions. Therefore, 
it is still needed to make a comprehensive study on the age distribution of star clusters 
over the entire disk of M82 to investigate its star cluster formation history. Since numerous 
previous studies on the same issue have led to various conclusions regarding the epochs of 
the encounter between M81 and M82, a study on the ages of the star clusters in the entire 
M82 disk, not restricted to some specific regions, will be of great merit. 

Another important issue is to find old halo clusters in M82 and to investigate their prop- 
erties. To date only a few old globular cluste rs were reported from a coup l e of spectroscopic 
observations of limited number of candidates (IKonstantopoulos et al.ll2009l : ISaito et al.ll2005f ) 
so that little is known about the old population in M82. Because of its proximity, M82 is 
an ideal target to search for old halo populations, if any. M82 is a very useful galaxy for 
understanding the existence and properties of old stellar populations in starburst galaxies. 

In this study, we present a survey of the star clusters in M82 using high resolution optical 
and near-infrared (U to H) images taken with H ST / ACS Wide Field Channel (WFC) and 
H ST /Wide Field Camera 3 (WFC3), covering not only the disk but also some halo area in 
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M82. The addition of U and H band data improves the capability to deal with the reddening 
and to reduce the age estimation uncertainty for star clusters. This paper is organized as 
follows: The methods of data reduction and analysis are described in Section 2. §3 presents 
the main results. Implication of these results are discussed in §4, and the primary results 
are summarized in §5. 

We adopt a distance to M82, 3.55 Mpc ((m — M) = 27.75 ± 0.07) determined from 
the me asurement of the /-band magnitude of the tip of the red giant branch (TRGB) mag- 
nitude Lee et all 120131 ) . The distance to M82 is similar to that t o M81, 3.63 ± 0.14 Mpc 



((m - M) = 27.80 ± 0.08) flDurrell. Saraiedini. fc Rupali I l2010l ; iGerke et al. 1 120111 ). At 



this distance, one arcsecond (arcmin) corresponds to 17 pc (1.02 kpc) in physical scale. The 
foreground reddening toward M82 is E(B — V) = 0.159 and the corresponding extinctions 
ar e A u = 0.862, A R = 0.685. A v = 0.526, A x = 0.308, A YJ = 0.177, and A H = 0.091 



( jSchlegel. Finkbeiner. fc David 119981 ) . The total magnitudes, colors, and position angle of 
M82 are B T = 9.3 ± 0.09, (B T - V T ) = 0.89 ± 0.01, and 65 degrees Jde Vaucouleurs et al. 



199ll ). At the adopted distance (without the internal extinction correction), the abso- 
lute magnitudes are Mg = —19.13 and My = —19.87. T he inclination angle of M82 
is i = 81.5 degrees in th e V band (ILvnds fc Sandaed Il963l ). i = 73 and 77 ± 3 degrees 



in the H and Ks band (Ichikawa et al. 



1! 



19951 ). and 77 ± 3 degrees in the JHKs images 



( IMayya. Carrasco. fc Lunall2005f ). respectively. We adopt 77 ± 3 degrees for analysis in this 
study. We use the center position of M82 from NASA/IPAC Extragalactic Database based 
on the Jordell-VLA Astrometric Survey (JVAS) and Cosmic Lens All-Sky Survey (CLASS). 



2. Data Reduction & Analysis 



2.1. Data 



Table [TJlists a summary of the data used in this study. We used the image data obtained 
by the Hubble Heritage Team using the H ST /ACS in F435W, F555W, and F8UW filter 
bands (called B, V, and I hereafter) (proposal ID 10776, PI: M. Mountain). In each filter 
band, six exposures were taken in a 3x2 mosaic configuration with the accumulated exposure 
times of 1600, 1360, and 1360 seconds, respectively. The final field coverage is about 12' x 8' 
centered on M82, covering the optical disk of M82. All the necessary data reduction processes 
were done by the Hubble Heritage Team, including the ACS image data reduction with the 
standard ACS pipeline (CALACS) and the image combination with MultiDrizzle task . Mor e 
detailed information on the observation and data reduction is given in lMutchler et al.l ( 120071 ). 
The full width at half maximum (FWHM) of a point source in the \^-band image is about 
0.1 arcsec. 
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We also used the ultraviolet (UV) and near-infrared (NIR) image data taken with 
H ST /WFC3 (Proposal ID 11360, PI: R. O'Connell) that partially covered the M82 disk 
in F336W (U), F110W (YJ), and F160W (H) bands. The full width at half maximum 
(FWHM) of a point source in NIR images is about 0.25 arcsec. All the image data provided 
from the HST archive are processed through the 'on-the-fly' calibration based on the best 
reference frames available. We have combined three exposures for U and two exposures for 
YJ and H using SWARP (jBertin et al.ll2002l ) to make a mosaic image in each filter band. 
The final mosaic images cover about 5' x 2'. See Figure [T] for the areal coverage of each 
band. 



2.2. Source Detection and Photometry 

Since M82 is an almost edge-on disk galaxy, it is very challenging to detect sources in 
the disk and central region of the galaxy. We combined the B, V, and I band images, and 
applied a ring median filtering method (with filter radii r in = 15 and r out = 20 pixels) with 
RMEDIAN task in IRAF to the combined image to generate a map of smoothly varying 
background brightness. Then we subtracted it from the original combined image, obtaining 
a master detection image. This method im proves the dete c tion e fficiency for the sources in 



the disk. A similar method was adopted by iHwang fc Led (120081 ) success fully to detect the 



star c lusters buried in the spiral arms of M51. We used Source Extractor (jBertin fc Arnouts 



19961 ) on the master detection image to detect sources with a threshold of 3cr and a minimum 



contiguous detected area of five pixels. 

We used PHOT task in IRAF/DAOPHOT package to derive photometry in UBVI, YJ, 
and H band images at the coordinates of detected sources. A single aperture with fixed size 
of r = 0.35" was adopted in all bands for the derivation of consistent colors. The aperture 
size was determined based on the FWHM of point sources in the NIR WFC3 image data 
with FWHMfs 0.3", which provides the poorest resolution in our data set. The photometric 
zero points for HST / ACS and WFC3 in Vega magnitude system are adopted from the STScI 



websites (http: //www. stsci.edu/hst/acs/ analysis/zeropoints/#tablestart 



and http://www.stsci.edu/hst/wfc3/phot_zp_lbn, respectively). 

The aperture magnitude requires correction for the lost light fraction due to the fi- 
nite size of an aperture. The aperture correction for point sources is determined based on 
the instrumental characteristics of the HST. However, for extended sources, the aperture 
correction also depends on their light profile characteristics and sizes. Therefore, we use 
an aperture c o rrecti on method for extended sources with a circular shape, as adopted by 



Hwang fc Led ( 120081 ). based on the correlation between the slope of light profiles and the 
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calculated value for the aperture correction. This correlation is obtained by calculating the 
slope of the aperture magnitude profiles between r = 5 pixels (= 0.25") and r = 7 pixels 
(= 0.35") for well-isolated and extended sources with circular shapes. The aperture cor- 
rection values are calculated by difference between aperture magnitudes with r = 7 pixels 
(= 0.35") and r = 15 pixels (= 0.75"). There is a good correlation between the slope and the 
aperture correction value (AMag= 2.94xslope— 0.13, rms= 0.09). To avoid over correction, 
we restrict the maximum correction value as 1.5 mag that is the largest aperture correction 
value of well-isolated sources. These aperture correction values are estimated only on the 
l/-band image to derive the l/-band total magnitudes of the star clusters. 



2.3. Star Cluster Selection 



Since selecting star clusters in nearby galaxies is a challenging process even with the 
high reso lution HST images, variou s methods have been used for the HST images of nearby 
galaxies. IScheepmaker et al.l (120071 ) selected star clus ters in M51 based on the source sizes. 
They measured the size of the sources using ISHAPE fjLarsenlll999[ ). and selected star cluster 
candidates with eff ective rad ii r e ff > 0.5 pc. Then they manually inspected the candidates for 
final selection. iMayya et al.l (120081 ) selected star clusters in M82 using the parameters pro- 
vided by Source Extractor. They checked FWHM, isophotal area, ellipticity, and magnitude 
of the detected sources to select star cl uster candidate s . The y also finalized the star cluster 
candidates through visual inspection. iHwang fc Led ( 120081 ) selected star clusters in M51 
using a combination of automatic screening and visual inspection. They used the FWHM 
and ellipticity information of the sources to define the cluster candidates. Star clusters were 
selected by inspecting radial profiles and 2D contours of the candidates. Both the automatic 
method and the visual inspection method are successful in cluster finding. However, the au- 
tomated methods have difficulties in finding the clusters located in the varying background. 
It needs additional visual check. Star clusters in M82 are often located in the bright and 
crowded background (especially in the disk), so that visual inspection is an indispensable 
step to find star clusters in M82. Therefore we use a combination of automatic screening 
and visual inspection to search for star clusters in M82 in this study. First we made a list 
of about 14000 candidates for visual inspection, selecting the sources with FWHM> 2.5 pix 
and V < 23 mag (after aperture correction). Then we visually checked the images, radial 
profiles, and isophote contour maps of the sources. Finally we selected 1105 star clusters, 
satisfying the following criteria: (1) round appearances, (2) extended radial profiles, and (3) 
lack of smooth extended structures. 



We compared our photometry of the M82 star clusters with those in previous stud- 
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ies ( jde Grijs. Bastian. fc Lamersl 120031 ; iMayya et al.l 120081 ; iKonstantopoulos et al.l 120091 ). as 
shown in Figure. [2j The (B — V) colors in these four studies show a ver y good agreement. V 
magnitudes in this study show an approximate agreement with those in|Mawa et al. fl2008h . 
but are systematically brighter t han those in other studies (jde Grijs. Bastian. &: Lamers 



2003 



Konstantopqulos et al.ll2009l) . This is due to the fact that the magnitudes in this study 



and IMayya et al.l (120081 ) provide total magnitudes for star clusters, while the other studies 
provide magnitudes for the fixed apertures. The mean differences derived with 2 sigma clip- 
ping are AV" = -0.737 ±0.049, A(B - V) = 0.025 ± 0.012, and A(V - 1) = 0009 ± 0.016 for 



de Griis. Bastian. fc Lamersl (120031 ) minus this st udy, A^ = -0.102 ± 0.018 , A(B - V) = 
0.002 ± 0.003, and A(V - I) = 0.035 ± 0.006 for iMavva et~aD ( 120081 ) minus this study, and 
A^ = -0.474 ± 0.043 , A(B - V) = -0.001 ± 0.007, and A ( V - I ) = -0.008 ± 0.016 for 



Konstantopoulos et al.l (120091 ) minus this study. IMayya et al.l ( 120081 ) adopted the isophotal 



magnitudes provided by Source Extractor, which are similar to the total magnitudes. 



2.4. Size Estimation of the Star Clusters 



To measure the size of star clusters, we us ed ISHAPE t hat was designed to estimate 
the size of slightly extended sources in images ( ILarsenl 119991 ) . It assumes that the image 
of a source can be modeled as convolution of a point spread function (PSF) and analytic 
functions for an extended source. We made a PSF image using DAOPHOT/SEEPSF task 
in IRAF for selected bright isolated stars. We adopted a MOFFAT15 model for analytic 
function. ISHAPE provides FWHMs of the sources, and we converted them to effective radii 
by multiplying 1.13, according to the guide of the ISHAPE manual. 



2.5. Estimation of Age, Extinction, and Mass of the Star Clusters 



We have estimated ages, extinctions, and masses of the star clusters in M82 by compar- 
ing the photometric data with the model spectral energy distributions (S EDs). Model SEDs 
are de rived from the simple stellar population (SSP) synthesis models by lBruzual fc Chariot 
( 120031 ). The models span ages from 1 Myr to 15 Gyr for metallicity range of Z= .0001—0.05. 



We assume a Salpeter initial mass function (IMF) of N(m)dm oc m~ 2 35 (ISalpeterlll955h with 
a lower mass limit with 0.1 M & and an upper mass limit with 100 M & . We adopted a solar 
metallicity for the star clusters in the disk based on the abundance of HII regions in M82 
( ISmith et al.l 120061 ). Note that using a Salpeter IMF could overestimate the cluster mass 
by a factor of 1.70-3.46 although r elative mass distribution of star clusters is not changed 
( jde Grijs. Bastian. &: Lamersl 120031 ). We adopted the Salpeter IMF, because most of the 
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previous studies adopted the same IMF and because we wanted to compare our results with 
them. We checked the effect of adopting the Chabrier IMF in comparison with the Salpeter 
IMF for SSP models in the (U — B) — (B — V) color-color diagram. The colors of the SSP 
models with these two IMFs are almost same for log t > 8, and they change by 0.1 for log 
t < 7. Therefore our results are affected little by adopting different IMFs. 

SED fitting based on broad band photom etry is known to hav e some limits, but is still 



very useful for estimating ages of star clusters. I Anders et al.l ( 12004 ) showed that at least four 



band data are needed to obtain reliable age estimation with an SED fit method. Especially, 
U and B bands are important for age estimation of yo ung star clusters , becau se those bands 



cover the Balmer jump that is a sensitive age indicator. iMaiz Apellanizl (120091 ) also suggested 
that the combination of U, optical, and NIR band data is useful to determine the ages of 
star clusters when the extinction of each cluster is unknown. M82 is almost edge-on, and 
most star clusters in the disk suffer from significant extinction. We utilize six UBVIiY J)H 
photometric bands for age and extinction estimation of the star clusters in the disk region. 
About 870 star cluster candidates in the central region of M82 were covered in six bands. 
However, only BVI images are available for about 200 star clusters, which are located mostly 
in the outer region of M82. 

The estimated ages of the clusters are not affected by the cluster size estimation, because 
the cluster colors were derived from the fixed aperture photometry. It is also noted that the 
estimated masses of the clusters do not depend on the cluster size estimation, because they 
were derived from the l/-band total magnitudes. 



3. Results 

3.1. A Catalog of the Star Clusters in M82 

In Table[2]we list a catalog of 1105 star cluster candidates (hereafter star clusters) with 
V < 23 mag. Brief descriptions on the data columns of the catalog are given as follows: 
Column 1 lists an identification number of a star cluster; Columns 2 and 3 list J2000.0 right 
ascension (RA) and declination (Dec) in degrees; Columns 4 lists aperture corrected \^-band 
magnitudes and corresponding errors; Columns from 5 to 9 list (B — V), (V — I), (U — B), 
(I — YJ), and (YJ — H) colors and corresponding errors for an aperture radius r = 0".35; 
and Column 10 and 11 list effective radii and ages, respectively. 
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3.2. Spatial Distribution of the Star Clusters 

Figure |3] displays the spatial distribution of the star clusters in M82. We defined the X 
and Y coordinates, which are projected distances along the major and minor axes, respec- 
tively. X and Y increase towards the northeast and the northwest, respectively, and are given 
by X = x' sin(PA) + y' cos(PA) and Y = — x' cos(PA) +y' sin(PA) where x' = (a — a ) cos(<5), 
y' = 5 — So, «o and 5o are the central coordinates of M82, and PA is the position angle for 
the major axis. 

Several features are noted in this figure. First, the most distinguishable feature is an 
elongated structure, showing that most star clusters are distributed on the disk: they are 
mostly located within V from the major axis and 3' from the minor axis. The distribution of 
the disk star clusters appears to show a discontinuity at ~ 3' along the major axis, indicating 
the disk star cluster populations may include two kinds. So we divided the disk area into 
two: the disk region (-3' < X < 3', -52" < Y < 52") and the disk-edge region (\X\ > 3', 
—52" < Y < 52"), as marked in Figure Second, there are a small number of star clusters 
scattered far from the disk. These are probably old globular clusters belonging to the halo 
of M82. Third, in the nuclear region that is known to have a re cent starburst , we fi nd 132 



star clusters inside the circular area with a radius 25" (425 pc). iMayya et al.l (120081 ) found 
260 star clusters in the same region, including fainter clusters in K-band images as well as 
clusters in very complex regions. Fourth, we find 80 star clusters in a 35" x 40" (595 pc x 
680 pc) field centered on the region B, a famous post-starburst region. We divide the entire 
field into four main regions for analysis: the main disk region, the nuclear region, the disk 
edge region, and the halo region as marked in Figure [3j Additionally, we select a sub region 
which is called region B for comparison with previous studies. This region is included in the 
disk region. 

We derived a number density of star clusters as a function of distance from the major 
axis (Y) to investigate the vertical structure of the star cluster system, as plotted in Figure 
HI The number density distribution appears to have a Gaussian shape with asymmetry. The 
peak is shifted from the center to the northern-side (Y > 0), and the number density of 
southern-side (Y < 0) is smaller than that of northern-side. This offset and asymmetry is 
considered to be due to the extinction and incompleteness of cluster detection. This shows 
that the star clusters in the southern-side suffer from more extinction than those in the 
northern-side, and that the southern-side (Y < 0) is closer to us than the northern-side. 

We fit the data with a hyperbolic s ecant function, iV = Nosech 2 (^-) where h z is a scale 



height ( Ide Grijs fc van der Kruitl Il996l ). plotting the result by the solid line in the figure. 
The fitting is performed for the range of Y > 0.2' to avoid the incompleteness problem. The 
vertical distribution of the star clusters is fit very well by the hyperbolic secant function for 
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0.2' < Y < V . We also investigate the vertical structure of stellar populations using the / 
band image. We sum counts with a range of —50" < X < 50" along the X-axis to derive a 
stellar light distribution along the Y-axis. This distribution of stellar light shows asymmetry 
similar to that of star clusters. We estimated the scale height of the vertical distribution 
of stellar light by fitting with a double hyperbolic secant function (for thin and thick disks) 
for the range of 10" < Y < 100". Derived scale heights of the thin and thick disks are 
h z (thin, star) = 8".38 ± 0".03 (-143 ±0.51 pc) and h z (thick, star) = 17".29 ± 0".05 (-294 
±0.85 pc), respectively. Thus the thick disk is twice thicker than the thin disk. 

We derived a value for h z of the star cluster distribution, 9". 64 ± 0".40 (164 ± 7 pc). 
This value is a little larger than that of the thin disk derived from stellar light, but this 
distribution includes the intermediate-age and old star clusters that are regarded as thick 
disk or halo populations. The genuine h z of the star clusters in the disk could be similar 
to that of stellar thin disk. It suggests that most star clusters are distributed similarly to 
stellar populations in the thin disk. 

It is also noted that there are several star clusters at Y > V above the fit function. 
These clusters may be a separate population from the disk population. They are considered 
to belong to the halo of M82. We defined the halo region by \Y\ > 3/i 2 (thick) (= 51". 87 = 
882 pc) derived from the /-band stellar luminosity profile for analysis, as marked in Figure 

El 

3.3. Photometric Properties of Star Clusters 

3.3.1. Color- Magnitude Diagrams and Color Distribution of the Star Clusters 

Figure [5] displays V — (B — V) color-magnitude diagrams (CMDs) of the star clusters 
in M82, and Figure [6] shows corresponding (B — V)q color histograms (for all star clusters 
and only bright star clusters with Vo < 21.5 mag): (a) the entire area, (b) the disk region, (c) 
the disk-edge region, (d) the halo region, (e) the nuclear region, and (f) the region B. Here 
magnitudes and colors represent the values dereddened according to the foreground redden- 
ing. Distinguishable features are as follows. First, the star clusters in M82 are dominated 
by the star clusters in the disk. They are composed of a dominant blue population with a 
color peak at (B — V)o ~ 0.45, and a weaker redder population that are mostly redder than 
(B — V)o ~ 1.0. We divide the cluster sample according to their color range: (B — V)o S 0.8 
(blue clusters), 0.8 < (B — V) ^ 1.3 (intermediate color clusters), and (B — V) > 1.3 (red 
clusters) for further analyses. Second, the CMDs and color distributions vary depending 
on the region. Third, the star clusters in the disk-edge region and the region B are mostly 



- 11 - 



blue. Fourth, the brightest star clusters in the disk-edge region, Vq ~ 20.5 mag, are much 
fainter than those in other regions. Fifth, the star clusters in the nuclear region include a 
large number of brighter star clusters compared with other regions, and they show a large 
scatter in color, indicating a presence of high differential reddening. Finally, the star clusters 
in the halo region are mostly red with 0.4 < (B — V)o < 1.2. These may be old globular 
clusters, showing the existence of an old halo in M82. A few of them have very red colors 
(B — V)q > 1.4, indicating that some of them may be reddened star clusters or background 
galaxies. 



3.3.2. Spatial Distribution of the Star Clusters with Different Colors 

The spatial distributions of the star clusters with different color ranges are shown in 
Figure [7J The blue star clusters ((B — V)o < 0.8) are mostly distributed in the disk, and 
several of them are found in the halo region. There are 69 and 48 blue star clusters in the 
region B and nuclear region, respectively. The intermediate-color clusters (0.8 ^ {B — V)$ < 
1.3) are also found mostly in the disk region and a few are in the halo region. Interestingly 
the distribution in the left is much wider than that in the right. It could be a result due to 
inhomogeneous extinction. A more detailed description is given in Section 3.6. Ten and 52 
intermediate-color clusters are found in the region B and nuclear region, respectively. The 
red star clusters {{B — V) ^ 1.3) are mostly located in the disk, especially in the right side of 
the disk. Only one red cluster is found in the region B, while 32 red clusters are detected in 
the nuclear region. Interestingly the red clusters are located along a narrow linear structure. 
We checked the Y J and H band images, and found the existence of dust lanes in the nuclear 
region. These dust lanes obscure the star clusters in the nuclear region, and the star clusters 
in the edge of the dust lanes are only shown as a narrow linear structure. The star clusters 
in the region B are usually blue, while more than half of star clusters in the nuclear region 
are red. It could be caused by the different age of star clusters or differential reddening. 
The region B is regarded as a window with less extinction and the nuclear region is usually 
considered as a high extinction region, so that the red colors of most clusters in the nuclear 
region may be a result of internal reddening. The linear distribution of the red star clusters 
in the nuclear region implies that there is a narrow window in the nuclear region or the 
clusters are formed in such a location. 
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3.3.3. Color- Color Diagrams of the Star Clusters 

Figure [8] displays the (U — B) — (B — V)o color-color diagrams (CCDs) of the bright 
star clusters with Vq < 21.5 mag in M82: (a) the entire area covered by U images, (b) the 
disk region, (c) the halo region, (d) the nuclear region, and (e) the region B. Note that 
the halo region is partially covered by U images so that only six objects are plotted. We 



also plotted the SSP models with solar metallicity from the Padova group ( jMarigo et al. 



20081 ). Several features are noted in Figure |HJ First, the disk star clusters are mostly in the 
region for 8 < log(age) < 10 and E(B — V) ~ 0.3. Second, the halo star clusters are old 
(log(age)> 9). Third, the star clusters in the nuclear region are mostly 8 < log(age) < 9 
and 0.1 < EiB — V) < 1.0. However, there are some clusters that are severely reddened 
(0.2 < E(B — V) < 2.0 or more). Fourth, the region B star clusters are mostly located in 
the region for 8 < log(age) < 9 and E(B — V) « 0.2. 



3.4. Luminosity Functions of the Star Clusters 



The luminosity function of the star clusters in M82 is shown in Figure IHJ The luminosity 
function of the star clusters with — 10 < My < —8 in the disk is fit well by a power-law 



star clusters in other late-type galaxies (e.g. 


Whitmore et al. 


1999 




Hwang & Lee 


2008; 


Santiago-Cortes. Mavva. & Rosa-Gonzalez 


2010 


) 


Annibali et al 


( 


2011 


) suggested that the 



slope of the luminosity function of the star clusters in the spiral galaxies becomes steeper 
when the luminosity of the fitting range is increased. The slope of the luminosity function 
for the star clusters in the disk region in M82 also follows this trend. The faint end in 
the luminosity function for the star clusters in the disk region becomes flatter at Vq ~ 21 
mag, which is a result due to incompleteness. This flattenin g is also seen i n the luminosity 
function of M82 star clusters derived from the same data by iMayya et al.l ( 120081 ). 



The luminosity function of the star clusters with —8 < My < —6 in the disk-edge region 
is also fit by a power-law function with a = — 2.77 ± 0.21, which is a little steeper than that 
of the star clusters in the disk region. The brightest star cluster in the disk-edge region is 
much fainter than those in other regions, as mentioned at Section 3.3.1. This result suggests 
two conflicting interpretations that the bright clusters could not be formed in this region 
or bright clusters were disrupted in this region. The former is more probable, considering 
that more massive clusters form in denser environments. The luminosity function of the 
star clusters in the nuclear region is fit by a power-law function with a shallower slope of 
a = —1.85 ± 0.06. However, it is noted that there is an excess in the bright end of the 
luminosity function for the nuclear region. It suggests that the starburst phenomenon in the 
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nuclear region could form very luminous star clusters, but we should deal carefully with this 
result because of extinction and incompleteness. The luminosity function of the star clusters 
in the halo region cannot be fit by a power-law function. It is noted that there are bright 
clusters with My < —8 in the halo region that are not found in the disk-edge region, while 
the numbers of the star clusters in the halo and disk-edge regions are similar. 



3.5. Size Distribution of the Star Clusters 



Figure [TU] displays the distributions of effective radii of the star clusters in M82: (a) 
all, (b) the blue star clusters, (c) the intermediate-color star clusters, and (d) the red star 
clusters. The star clusters have sizes of < r c g < 14 pc, but mostly smaller than 8 pc. The 
size distribution shows a peak at r e ff ~ 2 pc, and it smoothly decreases as the star cluster 
size increases. The median and mean values of the star cluster sizes are 2.99 pc and 4.02 
pc for all star clusters, 3.06 pc and 4.04 pc for the blue star clusters, 2.95 pc and 4.22 pc 
for the intermediate-color star clusters, and 2.14 pc and 3.33 pc for the red star clusters. 
The mean value of M82 st ar cluster sizes deri ved in this study is consistent with the peak 
value r e ff =3.8 pc given by iMayya et al.l (120081 ) . There are no significant differences of sizes 
between the blue and intermediate color star clusters (mean values are 4.04 pc and 4.22 pc, 
respectively). However, the sizes of the red clusters (a mean value of 3.33 pc) are a little 
smaller than the blue clusters (4.04 pc), but their standard deviations are too large (~ 5) to 
insist the difference. 

Figure [TlT a) displays the sizes of the star clusters versus the galactocentric distance 
(Rgc)- There is seen no significant trend of the sizes depending on Rgc- However it is noted 
that large star clusters with r e g > 10 pc are found mostly in the inner region of the disk 
at Rgc < l'-5. Figure fTTT b) and (c) displays the sizes of the star clusters versus l/-band 
magnitudes and (B — V) colors, respectively. Large star clusters with r e g > 10 pc are mostly 
fainter than V ~ 20 mag, while small star clusters with r e g < 10 pc have a wide range of 
magnitudes. The (B — V) color of the star clusters does not show any significant dependence 
on the size of the star cluster. 

The peak value for M82 star cl usters, 2 pc, is simi lar to that of the s tar clusters in 
other galaxies: e.g., 2.27 pc f or M51 (IHwang fc Led 120081) . ~ 2.3 pc for M83 (IBastian et al. 



20121 ). and ~ 1.5 pc for M31 (jVansevicius et al.ll2009i ). A relation betwe en cluster radii and 



galactocentric distance is found for the star clusters of M51 and M83 (jScheepmaker et al. 



20071 ; IHwang &: Ledl2008l ; IBastian et al.ll2012l ). The star clusters in the outer region of those 
galaxies have larger sizes than those in the inner region of those galaxies. However, the star 
clusters in M82 do not follow this trend. We also checked the relation between the sizes and 
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distances from the major axis, finding no significant trend. 

3.6. Age, Extinction, and Mass of the Star Clusters 

Figure fT2l displays the age distribution of ~630 star clusters in M82: (a) the entire region 
covered by the U-ba.nd images, (b) the disk region, (c) the halo region, (d) the nuclear region, 
and (e) region B. The ages of the star clusters are derived from SED fitting with two filter 
combinations, UBVI(YJ)H (left panels) and UBVI (right panels). Several distinguishable 
features are noted in Figure [T27 a). First, there is a dominant population of intermediate-age 
star clusters having ages from 100 Myr to 1 Gyr with a peak at ~ 500 Myr (log(age) = 8.7). 
This is seen in both UBV I{Y J)H fit results and UBVI fit results. There appears to be a 
separate weaker peak at ~ 125 Myr (log(age) = 8.1). This is better seen in the left panel 
than in the right panel. Second, there is a small number of young population with < 10 
Myr. Third, there is a sizable population of old star clusters with > 10 Gyr. 

Figure fl"2T b) to (e) show the variation of the age distribution of the star clusters depend- 
ing on their positions. First, the age distribution of the disk star clusters is similar to that for 
all star clusters because the entire cluster system in M82 is dominated by the disk clusters. 
Second, the number of the halo star clusters with age estimates is small, but they are mostly 
older than 1 Gyr. Third, the star clusters in the nuclear region show three populations: a 
major intermediate-age population with a peak value at ~ 400 Myr (log(age) = 8.6), a rela- 
tively small number of very young population with < 10 Myr as well as old population with 
> 10 Gyr. Fourth, the star clusters in the region B are mostly intermediate-age clusters. 

Figure H3] displays E{B — V) and masses versus ages of the star clusters in the halo, 
disk, and nuclear regions. E(B — V) values of the star clusters in the halo region are mostly 
small with E(B — V) < 0.3. E(B — V) values of the star cluster in the disk and nuclear 
regions are mostly distributed in the range of < E(B — V) < 1.0. The reddening value 
appears to increase as cluster ages decrease. In particular, for the young star clusters in the 
disk and nuclear regions with age < 10 Myr the reddening values are 0.8 < EiB — V) < 2.0. 
This indicates that young star clusters are usually located in dusty area suffering from high 
extinction, or that ages of these clusters might be under-estimated due to overestimation of 
extinction. 

The masses of the star clusters in M82 have a range with 1O 3 M < M < 1O 7 M , but 
mostly 10 4 M Q < M < 10 6 M Q . The lower boundary of the cluster masses is defined by the 
survey limit of this study. Most massive star clusters with M > 1O 6 M are located in the 
nuclear region. Several clusters with log(age) > 8 have very high masses with M ~ 1O 7 M , 
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but it should be treated with caution because over estimation of the extinctions could lead to 
high mass. Young star clusters with ages < 10 Myr have mostly low mass M <~ 1O 4 ' 5 M . 
Old star clusters are mostly massive with 1O 5 2 M < M < 1O 6 5 M . 

Figure HH displays the mass functions of the star clusters in the disk and nuclear region 
with different ages. There are several noticeable features as follows. (1) Upper parts of 
the mass functions of the young star clusters with log(age) < 9 in the disk and nuclear 
region mostly follow power-law functions (N(m)dm oc mP\ (2) The mass of the youngest 
star clusters with log(age) < 7 in both the disk and nuclear regions have similar slopes of 
power-law functions with fidisk = —1.53 ± 0.08 and (3 nuc = —1.54 ± 0.22 (fitting ranges of 
1O 3 ' 5 M < M < 1O 5 - 5 M for the star clusters in the disk region and 1O 4 -°M < M < 1O 6 -°M 
for the star clusters in the nuclear region). (3) The mass functions of the intermediate-age 
star clusters with 8 < log(age) < 9 in both regions have significantly different slopes of 
power-law functions. (4) The mass functions of the old star clusters in the disk do not follow 
a power-law function. 

It is known that the mass function of the young star clusters in most disk galaxies 
follows a power-law functions with (3 ~ — 2 (e.g. for M51 young star c lusters (log(age) < 7.0), 



P = -2.23±0.34 with 1O 3 6 M < M < 1O 41 M ( IHwang fc LeelbOlOf ): and for the young star 
clusters in Antennae ga laxies (log(age) < 7.0), f3 = -2. 14 ±0.04 with 1O 4 M < M < 1O 6 M 



( IWhitmore et al.ll2010f ). but the mass functions of the youngest star clusters in M82 have 
a shallower slope with ~ —1.5. It may be understood as a shallower initial cluster mass 
function (ICMF) of starburst galaxies, but the number of the young star clusters is too small 
to interpret this result as a top heavy ICMF. It is also noted that the mass function of the 
intermediate-age star clusters in the nuclear region is shallower with /3 = —1.46 ± 0.10 than 
that of typical young clusters, as shown in Figure IT4T c). However, this mass function can be 
interpreted as a result of star cluster disruption or incomplete survey. 

There are several extremely massive star clusters (M > 1O 6 ' 5 M ) in the nuclear region 
with 8 < log(age) < 9. They are more massive than typical globular clusters. The existence 
of the high mass star clusters could be interpreted as a result of a starburst phenomena 
of M82 at th e epoc h of the galaxy interaction. According to the spectroscopic study by 



Lancon et al.l ( 120081 ). two luminous star clusters have dynamical masses of 10 6 — 10 6 6 M ( 



05 

which is slightly less massive than our results suggest (~ 1O 7 M ), but they are still super 
massive. 

Figure [TB] and dH] display variations of E(B — V) and ages along the major axis X and the 
minor axis Y for the M82 star clusters. E(B — V) values are on average higher in the south- 
west side, indicating that the north-east side of the galaxy is farther than the south-west 
side. This result gives a hint for understanding the wider distribution of intermediate-color 
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star clusters (0.8 ^ (B — V)o < 1.3) in the north-east side. The star clusters with the largest 
values of E(B — V) are located mostly in the nuclear region. Mean values of E(B — V) are 
0.3-0.5 in the disk region, and ~ 0.8 in the nuclear region. On the other hand, mean cluster 
ages show little variation along the X or Y axes. 

We tested the effect of reddening errors to the derived ages and masses of the star 
clusters. The errors of derived E(B — V) in our SED fitting for M82 clusters range from zero 
to about 0.1 with a mean value of 0.04. A reddening error of 0.1 for E{B — V) leads to 0.1 
dex difference in log(age) and negligible difference in mass. 



4. Discussion 



4.1. Comparison with Previous Studies 



Ages and reddening values of star clusters in galaxies are often derived from the compar- 
ison of a combination of broad band photometry and/or spectr a with SSP models. In Figure 
IT71 we c ompare our e s timat e s for age and reddening with th ose in lde Grijs. Bastian. fc Lamers 
(ESI); ISmith et al.l (120061 ); iKonstantopoulos et al.l ( 120091 ) . Our ages and reddenings show a 
correlation with others's but with some scatters and systematic o ffsets. Our ages are on ay er- 
age slightly larger (~0.7 dex) than IKonstantopoulos et al. (2009)' s andlSmith et al.l ( 120061 ) 's. 
but slightly smaller (~0.12 dex) than Ide Grijs. Bastian. &: Lamersl (2003)'s. Our r edden- 
ing values a r e on a verage slightly smaller (~0.2) tha n IKonstantopoulos et al.l (120091) ? s an d 
Smith et al. (2006)'s, but slightly larg er (~0.2) than de Griis. Bastian. fc Lamersl ( 2003 )'s 



that are mostly zero values. Note that Konstantopoulos et al. ( 20091 ) derived ages from the 
spectra first, then estimated reddenings from BVI colors based on these ages. 

These differences c ould be explained by raising a few points. The first possibility is 
the different filter sets. Ide Grijs. Bastian. fc Lamersl (120031 ) used BVIJH filter set, while 
we used UBVI(YJ)H. It is known that the existence of U band data will lead to a better 
determination of a ges and reddenings with SED fit methods, especially for young clusters. 
Smith et al.l (120071 ) found, from the UBVI photomet ry of the star clusters in the regio n 
B, and that their ages are much younger than those of Ide Grijs. Bastian. &: Lamersl ( 120031 ). 
while their redde nings are much larger than the latte r. This is consistent with our findings as 
above. Note that Ide Grijs. Bastian. fc Lamersl ( 120031 ) obtained almost zero reddening values, 
which may lead to slightly older ages than those of ours (~0.2 dex) with estimated reddening 
of E(B-V) -0.3. 



The secon d possibility i s the d iffere nt SSP models. We used S SP m odels from lBruzual &: Chariot 



( 120031 ). while ISmith et al.l (120061 ) and IKonstantopoulos et al.l (120091 ) adopted the STAR- 
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BURST99 models from iLeitherer et all (119991) . This model differenc e could lead to the dif- 
ferent results of the age estimation. Recently iRodrfguez-Merino et all (120111 ) investigated the 
evolution history of M82 using multi- wavelength photometry of unresolved stellar light in the 
disk. They used 5 arcsec apertures to derive spatially resolved star formation history of the 
M82 disk with multi-wavelength data from UV to NIR. They co mpared observational data 
with SED models from the Padova Group (IBressan et al.lll993l ) with solar-metallicity and 
Salpeter IMF. This model is very similar to the model adopted in this study. They derived 
mean values for age and reddening of the apertures for the region B, log(age) = 8.60 ± 0.06 
and E(B — V) = 0.44 ± 0.05. These values for integrated stellar light are similar to the 
mean values for the star clusters in the region B derived in this study, log(age) ~ 8.7 and 
E(B — V) — 0.30. Therefore the model differ e nce c ould be a main source of age differe nces 
between our results and those of ISmith et al.l (120061 ) and iKonstantopoulos et al.l ( 120091 ). 



The third po ssibility is the diffe rent methods. We used SEP fittin g with broad band 
photometry, while ISmith et al.l (120061 ) and IKonstantopoulos et al.l (120091 ) derived cluster pa- 
rameters from spectroscopy. Both methods yield similar results in general, but with a large 
scatter. It is not easy to disentangle clearly each possibility among these. 



4.2. Halo Star Clusters in M82 



There is only a little information on the existence of old stellar population in the halo 
of M82, while there are many observational reports showing the existence of ionized gas , 



atomic gas, and molecular clouds in the halo region of M82 dTavlor. Walter, fc Yunl 12001 



Chynoweth et al.l l2008t IVeilleux. Rupke. fc Swatersl 120091 ; iRoussel et al.l 120101 ) . Some ob- 
servational stu dies reported the existence of intermediate-age stars outside the disk of M82 
( lDavidgell2008l ). Only a small number of globular clusters in the field of M82 were discovered 
by pr evious spectroscopic studies (two in ISaito et al.ll2005l and two in IKonstantopoulos et al. 
20091 ). Thes e globular clusters were fo und only in the south of M82 in the direction toward 
M81 so that iMayya fc Carrasco J2009h suspected that they may not belong to M82. 



In this study, we have discovered 35 star clusters in the halo region of M82. The color 
distribution of these halo clusters show a dominant blue peak at (B — V) ~ 0.65 and they are 
mostly old (> 10 Gyr). In addition they are located in a wide area surrounding M82. These 
suggest that a significant population of old globular clusters does exist in the halo of M82. 
However, there may be some contamination due to background galaxies. Measurements 
of radial velocities of these objects are needed to distinguish the halo members and the 
background galaxies. 
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Another noteworthy point for the halo clusters in M82 is that there are some clusters 
whose stars are partially resolved. With the resolving power of H S T / ACS, it is possible to 
resolve some bright stars in M82 (ILee et al.ll2013l : I Jang et al.ll2012[ ). If we find star clusters 
that are partially resolved into stars in M82, then they can be considered to be the members 
of M82. We made a smoothed image for a field including each star cluster, and subtracted 
it from the original image to reveal any feature of resolved stars in each cluster. Then we 
inspected the subtracted images, finding 13 halo star clusters showing features of resolved 
stars, as shown in FigurelTSl 

We transformed the photometry of the halo star clusters from the HST filter system into 
Johnson-Cousins filter system using the calibration information from ISirianni et all (J2005) 
to co mpare the halo star clusters with the globular clusters in the Milky Way Galaxy ((Harris 
19961 ). Figure[T§ldisplays a (B— V) Johns on — (V— -On.johnson ffi ^ ne na ^° s ^ ar crus t ers - SSP 
models for Z = 0.02,0.004 and 0.0001 faruzual fc Charlotl liooih are also plotted. The star 
clusters with and without resolved stars are plotted in different symbols (circles and squares, 
respectively). Most halo star clusters have color ranges with 0.5 < (B — V r )o,j hnson ^1-2 and 
0-8 ~ (V ~ -Oo,Johnson ~ 1-5- This color range is similar to that of the Milky Way globular 
clusters. The resolved star clusters and unresolved star clusters have similar color ranges 
except for a few extremely red unresolved clusters. Therefore most of the unresolved star 
clusters are also probably the globular clusters of M82. Three reddest objects may be highly 
reddened star clusters. 



Since the halo region is far from the disk plane (projected distance > 880 pc), we can 
assume that the internal extinction is negligible. The halo population in late-type galaxies 
is usually metal poor so that we assume a lower metallicity for the halo star clusters in 
M82. Based on these assumptions, the halo star clusters mostly concentrated on the region 
at (B — V)o ~ 0.65 and (V — J) ~ 1.0 in Figure 19 are consistent with 10 Gyr age with 
Z= 0.0001. Several halo clusters are located in the region extended to the redder direction 
along the SSP model. These red halo clusters could be candidates of reddened old metal-rich 
star clusters. Four of them are located relatively closer to the disk (< lkpc), so that they 
might have been reddened due to internal extinction. 



4.3. Comparison of M82 Halo Star Clusters with M81 Globular Clusters 

We compared the M82 halo star clusters with the M81 globular clusters in the Vo,j hnson — 
(B — l / )o,johnson diagrams and (B — V)o t j \ mson color distributions in Figurel20l It is relatively 
easy to select halo star clusters in M82, because it is an almost edge-on galaxy. In contrast, 
it is difficult to select halo globular clusters in M81, because M81 is almost face-on and many 
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globular clusters were found in the projected plane of the disk region of M81. 



Nantais et al.l ( 1201 ll ) found 214 globular cluster cand idates in M81 including 85 con- 



firmed by spectroscopic studies (INantais &: Huchral 120101 ) . and provided a catalog of the 



globular cluster candidates with Johnson- Cousins BVI photometry. The brightest glob- 
ular clusters in M81 are about one magnitude brighter than those in M82. The color 
ranges of the M82 halo star clusters and M81 globular clusters are similar, with 0.5 < 
(B — V) 0i johnson ^ 1-3, although four of the unresolved halo star clusters in M82 have redder 
colors than (B — V) rs 1.3. 

The color distributions of the M82 and M81 clusters are different. The M82 clusters 
show an asymmetric distribution, while the M81 clusters have a Gaussian-like distribution. 
The peak color for the M82 clusters, (B — K)o,Johnson ~ 0.65, is bluer than that for the M81 
clusters, at (B — ^)o,johnson ~ 0.95. This difference is seen both in the total sample and in 
the confirmed globular cluster samples (resolved star clusters in M82 and spectroscopically 
confirmed globular clusters in M81). If the internal extinction for the clusters is negligible 
and the mean ages are similar, the difference in the peak color could be explained by different 
metallicities: the M82 halo clusters are on average more metal poor than the M81 globular 
clusters. However it is noted that there are a small number of M82 halo clusters that have 
similar colors to the peak color of the M81 globular clusters. These results show that M82, 
showing a recent starburst, have an old stellar halo including metal poor old globular clusters, 
In addition, the origin of the metal-poor halo clusters in M82 appears to be different from 
that of the M81 globular clusters. 



4.4. Evolution of M82 



M82 is a famous archetype of the starburst galaxy, and many studies have concentrated 
on understanding of the recent starburst. Previous studies suggested that there might have 
been two starbursts in the central region of M82: one at ~1 Myr ago and another at ~5 
Myr ago (IForster Schreiber. Genzel. Lutz. fc Sternberg! 120031 ). The first starburst started in 
the nuclear region as well as in the bar, while the second burst in the circum-nuclear region. 



Recently, iMayya fc Carrascol (120091 ) noted from the analysis of stellar light that M82 
differs from normal late type galaxies: normal spiral galaxies mostly have disks with con- 
tinuing star formation, and its mean stellar age is old (> 5 Gyr), while the mean age of 
the disk of M82 has relatively young ages (~500 Myr) and star formation was quenched 
at about 100 Myr ago. This is consistent with the results derived from the star clusters in 
this study. Young clusters in M82 have been used to study the impact driven by dynamical 
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interactions with M81. IChandar. Tsvetanov. &: Fordl (120011 ) found 114 star clusters in eight 
WFPC2 fields of M81. They showed that the young cluster formation in M81 started at 
~600 Myr ago from the analysis of BVI CCD. This is consistent with the population of star 
clusters with an age peak at ~500 Myr in M82 derived in this study. The existence of almost 
coeval star clusters in M81 and M82 suggests that the interaction between two galaxies at 
this epoch may have enhanced the star cluster formation in both galaxies. 



Mayya fc Carrascol (120091 ) suggested three possible scenarios to explain the star forma- 



tion history derived from stellar light in the disk as follows: (1) a starburst galaxy formed by 
flyby encounter of low surface brightness (LSB) galaxies, (2) a rejuvenated disk galaxy with 
a truncated stellar disk, and (3) a newly formed galaxy. They supported the first scenario, 
but they could not rule out the other scenarios, either. We discuss each scenario in the 
following considering the results derived in this study. 



The first scenario is based on theoretical simulations (IMihosI Il999l ). This simulation 
suggested that flyby encounter of a companion galaxy could have induced the disk-wide star 
formation in the host gas-rich LSB. We found that a significant star formation over the entire 
disk at about 500 Myr ago from the age distribution of M82 disk clusters as shown in Figure 
If we adopt this model for the cause of the observed disk-wide star formatio n, then the 



epoch of the interaction between M81 and M82 should be about 500 Myr ago. lYun et al. 



( 119931 ) suggested that the gas around M82 may have originated from M82, and it means 
that M82 was a gas- rich galaxy before its interaction with M81. These results support this 
scenario. However, this scenario cannot explain the recent central starburst in M82. 



The second scenario was proposed by ISofud (119981 ). In this scenario, the disk of M82 
was stripped by the interaction and a new disk was formed. Then the old stars stripped 
from the disk and halo should be discovered in th e M81 group. These remnants of stripping 
were not yet found, but iMayya fc Carrascol (120091 ) noted that the stripped stars may be low- 
mass main sequence stars and they are fainter than the detection limit of HST observation. 
However, stel lar ages and velocity fields o f M82 are not consistent with the expectation of 
this scenario ( IKonstantopoulos et al.ll2009l ) . This scenario suggested that the stellar contents 
of M82 are dominated by the bulge population. However, we found that the star clusters 
with ages ranging from very young (< 10 Myr) to very old (> 10 Gyr) exist in the nuclear 
region. It suggests that the star clusters in the nuclear region are not dominated by the 
bulge population. 

The last scenario is that M82 is a newly formed galaxy. However, we found the existence 
of the old globular clu ster system in this study, and AGB stars are also found in the halo 
region (jDavidgel 120081 ). These show that M82 is not a newly formed galaxy. Therefore we 
can rule out the second and third scenarios. Our results support the first scenario, but the 
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existence of recent starbursts in the nuclear region remains to be explained. 

5. Summary 

In this study we presented a photometric study of the star clusters in the starburst 
galaxy M82 based on high resolution HST images. We found 1105 star clusters in the entire 
H ST /ACS field of M82, and derived their photometry from UBVI(YJ)H images. Primary 
results are summarized below. 

• A majority of the star clusters are located in the disk of M82. The scale height of their 
vertical number density profile is derived to be h z = 9". 64 ± 0".40 (164 ± 7 pc), which 
is similar to that of the thin stellar disk of M82. We also found 35 star clusters in the 
halo region of M82. 

• The CMDs and color distribution of the star clusters vary depending on the regions. 
The star clusters in the disk of M82 are composed of a dominant blue population 
and a weaker red population. The blue populations in the disk have a color peak at 
(B — V)o w 0.45, and the red populations in the disk are mostly redder than (B — V) « 
1.0. The star clusters in the disk-edge region are mostly bluer than (B — V) ~ 0.8, 
while those in the nuclear region have a large color range of 0.2 ;$ (B — V) ;$ 2.0. 
The star clusters in the halo region have a color range of 0.4 < (B — V)o J$ 1.5 that is 
similar to that of the globular clusters in the MWG. 

• The luminosity function of the disk star clusters follows a power-law function with 
a = —2.04 ± 0.03 that is similar to those of young clusters in other galaxies, while the 
luminosity function of the star clusters in the halo and nuclear region do not follow 
power-law functions. The luminosity function of the star clusters in the disk-edge region 
is also fit by a power-law function but with a slightly steep slope, a = —2.77 ± 0.21. 

• The size distribution of the star clusters has a strong peak at ~ 2 pc, and has a long 
tail to larger radii. There is no significant trend between the star cluster size and 
galactic central distance. 

• We estimated the ages of about 630 star clusters in M82 with U BVI(Y J)H band data, 
and found that the most dominant population shows a peak at ~ 500 Myr. There are 
also some star clusters with younger (< 10 Myr) and older (> 10 Gyr) ages. 

• The mass function of the young star clusters (log(age)< 7.0) is fit by a power law with 
a slope (/3 ~ —1.53 ± 0.08) shallower compared with that of other late type galaxies 
such as M51 and Antennae galaxies. 
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• We found two populations of star clusters in the nuclear region: very young clusters 
with ages of 5 Myr and intermediate age clusters with ages of ~500 Myr. The former 
represents the recent starburst, while the latter are consistent with those in the disk. 
The brightest young clusters have high masses (> 1O 5 M ), indicating that massive 
clusters form preferentially in the high density regions. 

• The colors and ages of the halo star clusters in M82 are similar to those of the globular 
clusters in the Milky Way Galaxy, so that they could be regarded as old globular 
clusters in M82. It suggests that M82 has an old stellar halo that is similar to that of 
other normal galaxies. From the comparison of (B — l / )o,j hnson colors, we have found 
that the M82 globular cluster system is on average more metal poor than M81 globular 
cluster system. 

• Our results suggest that there were four epochs of starburst in M82. The first epoch is 
when the old globular clusters were formed. The second is when a majority of the star 
clusters were formed in the entire disk, which started 1 Gyr ago and stopped 100 Myr 
ago. The third is when star clusters were formed in the nuclear region < 10 Myr ago. 
The fourth is when star clusters were formed recently in the central region (maybe ~3 
Myr ago). 



This work was supported by the National Research Foundation of Korea (NRF) grant 
funded by the Korea Government (MEST) (No. 2012R1A4A1028713). 
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Table 1: A summary of data. 



Filter 


Proposal ID 


Instrument 


Exposure times (s) 


F435W(B) 


10776 


HST/ ACS WFC 


1600 


F555W(V) 


10776 


HST/ ACS WFC 


1360 


FSUW(I) 


10776 


HST /ACS WFC 


1360 


F33QW(U) 


11360 


H ST /WFC3 UVIS 


1050, 1215, 1620 


FU0W(YJ) 


11360 


iJ,ST/WFC3 IR 


598, 1195 


FIQOW(H) 


11360 


HST /WFC3 IR 


598, 2395 



Table 2: A catalog of the star clusters in M82 



ID 


R. A. (,72000) 


Dcc.(J2000) 


V 


(B-V) 


(V - I) 


(U-B) 


(/ - YJ) 


(YJ - H) 


r cff [pc] 1 


og(agc) [yr] 


481 


9 55 47.07 


69 40 42.3 


16.289 ± 0.007 


1.127 ± 0.012 


1.531 ± 0.009 


0.764 ± 0.015 


1.050 ± 0.013 


0.697 ± 0.019 




8.67 


611 


9 55 52.68 


69 40 46.0 


17.061 ± 0.028 


1.032 ± 0.044 


0.906 ± 0.092 


-0.139 ± 0.038 


0.545 ± 0.353 


1.140 ± 0.368 




9.17 


628 


9 55 53.47 


69 40 51.1 


17.125 ± 0.004 


1.322 ± 0.006 


2.061 ± 0.013 


0.500 ± 0.009 


1.166 ± 0.051 


0.815 ± 0.069 


3.23 


8.25 


624 


9 55 53.34 


69 40 49.8 


17.207 ± 0.020 


1.128 ± 0.021 


1.425 ± 0.075 


0.382 ± 0.011 


1.162 ± 0.194 


0.738 ± 0.260 


4.76 


9.69 


561 


9 55 50.42 


69 40 47.3 


17.407 ± 0.013 


1.260 ± 0.017 


1.521 ± 0.047 


0.366 ± 0.016 


1.267 ± 0.098 


0.618 ± 0.139 


4.91 


9.38 



to 

00 



(This table is available in its entirety in a machine-readable form in the online journal. A portion is shown 
here for guidance regarding its form and content.) 
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Right Ascension 



Fig. 1. — A grey scale map of the H ST/ ACS /-band image of M82. Solid, dotted, and dashed 
lines represent the areas covered by HST/ ACS (FA35W(B), F555W(V) , and F8UW{I)), 
HST/WFC3 UVIS (F336W(U)), and NIR (FllO W(YJ) and F I GOW (H)) images, respec- 
tively. A large ellipse represents the survey area of iMayya et al.l ( 120081 ) with 250" of major 
axis, 60" of minor axis, and 65 degrees of position angle. Arrows at the bottom-left corner 
represent X and Y directions defined in Section 3.2. 
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16 18 20 22 16 18 20 22 16 18 20 22 

V (This study) 

Comparison of our photome t ry of M82 star clusters with those in previ ous stud 



de Griis. O'Connell. fc Gallagher! (120011) (upper panels), iMavva et al.l ( 120081 ) (middle 
panels), and iKonstantopoulos et al.l (120091 ) (lower panels). Left, middle, and right panels 
represent the differences in V, (B — V), and (V — I), respectively. A represents this study 
minus previous studies. Filled circles represent the data used for deriving average values 
(dashed lines) with 2-sigma clipping. 
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Fig. 3. — Spatial distribution of the star clusters in M82. The large solid-line polygons 
indicate the coverage of H ST /ACS (BVI), H ST /WFC3 UVIS (U), and H ST /WFC3 NIR 
(YJ, H) images. A large circle in the center represent the nuclear region with a radius of 
25" (425 pc), and a box left of the circle marks the region B (35" x 40", 595 pc x 680 pc). 
The horizontal dashed lines represent the boundary between the disk region and the halo 
region. The dotted lines indicate the boundary between the disk region and the disk-edge 
region. Crosses, diamonds, filled circles, and open squares represent the star clusters in the 
disk, nuclear region, halo, and disk-edge region, respectively. One arcmin corresponds to one 
kpc. 
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300 




Y [arcmin] 

Fig. 4. — The number density profile along the minor axis (Y) for the star clusters in M82 
(histogram), which is fit with a hyperbolic secant function (solid lines). The vertical dotted 
line represent a lower limit for the fitting range. For comparison, fitting results for stellar 
light distribution are also shown for the thin disk (dotted line), thick disk (dot-dashed line), 
and sum of the two disks (dashed line). 
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Fig. 5. — Color-magnitude diagrams of the star clusters in M82: (a) the entire region, (b) 
the disk region, (c) the disk-edge region, (d) the halo region, (e) the nuclear region, and (f) 
region B. The bars at the left side in (a) show mean (B — V) errors for magnitude bins. 
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Fig. 6. — (B — V) color histograms of all star clusters (solid line histogram) and bright star 
clusters with Vo < 21.5 mag (dashed line histograms) in M82: (a) the entire region, (b) 
the disk region, (c) the disk-edge region, (d) the halo region, (e) the nuclear region, and 
(f) region B. Vertical lines in the panel (a) represent the boundary of three groups: blue, 
intermediate-color, and red clusters. 
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Fig. 7. — Spatial distributions of the star clusters in M82 with different color ranges: (a) all, 
(b) blue clusters ((£> — V) < 0.8), (c) intermediate color clusters (0.8 ^ (B — V) < 1.3), 
and (d) red clusters (1.3 ^ (B — V)o). The large solid-line box indicates the coverage of 
H ST /ACS (BVI). The dashed lines represent the boundary between the disk and halo 
regions, and the dotted lines represent the boundary between the disk and disk-edge regions. 
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Fig. 8. — (U — B) - (B — V)o color-color diagrams of the bright star clusters with Vo < 21.5 
mag in M82: (a) the entire region, (b) the disk region, (c) the halo region, (d) the nuclear 
region, and (e) region B. The solid curved lines sh ow the simple stellar population model for 
Z = Z Q (except for the halo, Z = 0.001) given by lMarigo et al.l ( 120081 ). shifted according to 
E(B — V) = 0.0, 0.2, 0.3, 1.0, and 2.0. The arrows represent reddening vectors with Ay = 9. 
Error bars at the right-upper corner of all panels represent median values of photometric 
errors. 
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Fig. 9. — Luminosity functions of the star clusters in the disk (a), disk-edge (b), nuclear (c), 
and halo (d) regions of M82. The thick solid lines in (a), (b), and (c) show a power-law fit 
to the luminosity function of the star clusters over the 18.0 > Vq > 20, 20.0 > Vq > 22, and 
17.5 > Vq > 19.5, respectively (the thick dashed lines represent an extrapolation). 



-38 - 



150 

£ 100 
50 

30 
£ 20 
10 



- 




i 

■ ^o\ aii 
(aj All 

i 
i 

tZJ— Median-2.99 pc 

< — Mean-4.02 pc 
i 

\ : 

i 

: ^ 




< 


(D) (B-V) < U.5 

— Median=3.06 pc 
< Mean-4.04 pc 

i 

■ i ■ ■ ■ i ■ iT|T!Th"i ■ i-r^H ■ 






1 ■ 1 ■ ■ ■ 1 ■ ■ ■ 1 ■ ■ ■ 1 ■ iTTTH ■, 
1 

C ^ 1 U . O _ ^ D - V jg N l.J 

1 
1 

^JL_Median=2.95 pc 

< — Mean-4.22 pc 

i : 

j— i 

i 

''in ! 






— Median=2.14pc 
< — Mean=3.33 pc 

• ■ 1 



2 4 6 8 10 12 14 2 4 6 8 10 12 14 
r e ff [pc] r eff [pc] 



Fig. 10. — Size distributions of the star clusters in M82: (a) all, (b) blue clusters ((B — 
V)o < 0.8), (c) intermediate color clusters (0.8 ^ (B — V) < 1.3), and (d) red clusters 
(1.3 ^ (B — V)o). The vertical solid and dashed lines indicate the median and mean values 
of effective radii, respectively. 
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Fig. 11. — Size versus galactocentric radii, Rqc ( a ) ; Vq (b), and (B — V) (c) for the star 
clusters in M82. Small filled circles, large open circles, and open squares represent the star 
clusters with good-fit sizes (x 2 < 10) in the disk, halo, and disk-edge region, respectively. 
Triangles with error bars in (a) represent mean values of the star cluster sizes in each bin. 
The error bars of X axis direction and Y axis direction show radial bin sizes and standard 
deviations of mean sizes, respectively. 
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Fig. 12. — Age distributions of the star clusters in M82: (a) the entire region, (b) the disk 
region, (c) the halo region, (d) the nuclear region, and (e) region B. The solid and dashed 
lines indicate all star clusters and bright star clusters with Vq < 21.5 mag, respectively. The 
left and right panels show the ages of the star cluster derived from UBVI(YJ)H and UBVI 
, respectively. 
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Fig. 13. — E(B — V) versus log(age) (a) and log M/M versus log(age) (b). Crosses, 
filled circles, and diamonds represent the star clusters in the disk, halo, and nuclear regions, 
respectively. 
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Fig. 14. — Mass functions of the star clusters in the nuclear region (solid lines) and disk 
(dashed lines) with different ages: (a) log(age)< 7, (b) 7 < log(age) ^ 8, (c) 8 < log(age) 
< 9, and (d) 9 < log(age). Linear lines represent power law fits for the upper mass range. 
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Fig. 15.— E(B — V) (upper panel) and ages (lower panel) versus X for the M82 star 
clusters. The squares with error bar represent the mean values of ages and E(B — V) for 
the star clusters with log(age)< 10. Data with log(age)=10.3 represent the star clusters for 
which fitting failed and that were assigned very old ages. 
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Fig. 16.— E[B — V) (upper panel) and ages (lower panel) versus Y for the M82 star 
clusters. The squares with error bar represent the mean values of ages and E(B — V) of the 
star clusters with log(age)< 10. 
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Fig. 17. — Comparison of ages (left panels) and E(B — V) (right panels ) deriv ed in this 
study with previous stud ies ((upper panels) Ide Griis. Bastian. &: Lamersl (120031 ) . (middle 
panels) ISmith et al.l ( 120061 ) . and (lower panels) iKonstantopoulos et al.l ( 120091 )). Dashed lines 
represent one-to-one relations. 
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Fig. 18. — Thumbnail J-band images of the halo star clusters that are partially resolved into 
individual stars. Each image covers 4" x 4". Smoothed images were subtracted from the 
original images for showing better resolved stars. 
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Fig. 19. — (B — V)o, Johnson — (V — I)o,johnson diagram of the star clusters in the halo region 
of M82. Here (B — V)o, Johnson and (V — /)o, Johnson are given in the Johnson-Cousins system, 
and foreground reddening is corrected.. The size of crosses represents an error of each color. 
Open circles and open squares indicate the resolved and unresolved star clusters, resp ectively. 
Filled circles represent the globular clusters in the Milky Way Galaxy (lHarrislll996f ). Dotted 
line, dashed line, and solid lines show the simple stellar p opulation models with Z = 0.0001, 
0.004, and 0.02, respectively, given in (IMarigo et al.ll2008l ). Numbers labeled along the model 
lines represent log(age). The arrow represents a reddening vector with Ay = 1. 
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(c) SC in M82 halo 
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Fig. 20. — Vo,j hnson — {B — K)o,johnson diagrams (upper panels) and color histograms (lower 
panels) of the st ar clusters in the ha lo of M82 (left panels) in comparison with the globular 
clusters in M81 (INantais et al.ll201l[ ) (right panels). Here Vo 5 j hnson and (V — /)o,johnson ar e 
given in the Johnson- Cousins system, and foreground reddening is corrected. Dots and 
asterisks in (a) represent unresolved and resolved star clusters in the M82 halo, respectively. 
Asterisks and dots in (b) indicate M81 globular clusters confirmed by spectroscopic studies, 
and good globular cluster candidates, respectively. Solid and dashed lines in (c) represent 
all star clusters and resolved star clusters in the M82 halo. Solid and dashed lines in (d) 
represent all globular clusters and spectroscopically confirmed globular clusters in M81. 



